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Abstract. Because of the tremendous densities that exist in the cores of neutron stars, a significant fraction 
of the matter in the cores of such stars is likely to exist in the form of hyperons. Depending on spin frequency, 
the hyperon content changes dramatically in rotating neutron stars, as discussed in this paper. 
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Rotating neutron stars are called pulsars. Three distinct 
classes of pulsars are currently known. These are (1) rota- 
tionpowered pulsars, where the loss of rotational energy 
of the star powers the emitted electromagnetic radiation, 
(2) accretion-powered (X-ray) pulsars, where the gravita- 
tional potential energy of the matter accreted from a low- 
mass companion is the energy source, and (3) magnetars, 
where the decay of a ultra-strong magnetic field powers 
the radiation. The matter in the cores of rotating neutron 
stars is compressed to ultra-high densities that may be 
more than an order of magnitude greater than the den- 
sity of atomic nuclei. This makes (rotating) neutron stars 
superb astrophysical laboratories for a wide range of fas- 
cinating physical studies . This includes the physics 
of hyperons in cold ultra-dense matter, whose thresholds, 
according to model calculations, are easily reached in the 
cores of neutron stars (see Fig. [T]), depending on the mass 
and rotational frequency of a neutron star. The most rapidly 
rotating, currently known neutron star is pulsar PSR J1748- 
2446ad, which rotates at a period of 1.39 ms (which cor- 
responds to a rotational frequency of 719 Hz) [4j. It is fol- 
lowed by PSRs B1937-h21 [5| and B1957-h20 whose ro- 
tational periods are 1.58 ms (633 Hz) and 1.61 ms (621 Hz), 
respectively. Finally, we mention the recent discovery of 
X-ray burst oscillations from the neutron star X-ray tran- 
sient XTE J 1739-285 [T], which could suggest that XTE 
J1739-285 contains the most rapidly rotating neutron star 
yet discovered. Rapid rotation changes the structure and 
composition of neutron stars dramatically, and leads to 
novel phenomena like frame dragging (Lense Thirring ef- 
fect), as discussed in this paper. 
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Fig. 1. Change of central density with rotational neutron star 
frequency 0. to = 140 MeV/fm'^ denotes the density of nuclear 
matter, /2k is the Kepler frequency, and M{0) is the star's 
mass at zero rotation. 



2 Stellar structure equations 

Since neutron stars are objects of highly compressed mat- 
ter, the geometry of spacetime is changed dramatically 
from flat space by these objects. Neutron star models are 
thus to be computed from Einstein's field equations of 
general relativity (/x, i/=0,l,2,3), 
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G^, = i?^, -i(/^,i? = 8^r^,(e,P(e)), (1) 

which couples Einstein's curvature tensor, G^i,, to the 
energy-momentum density tensor, T^,^, of the stellar mat- 
ter. The quantities R^, = r;<,_, - F^,,^ + F^^F^^ - 
r^g.r^j,, g^y, and R = Rp^i/g^" in denote the Ricci ten- 
sor, metric tensor and scalar curvature, respectively [2]. 
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(Commas followed by a Greek letter denote derivatives 
with respect to space-time coordinates, e.g. ,u = d/dx'^ .) 
The ChristofFel symbols are defined as = (1/2)5'^^ 
(5/iA, V + gv\, — g^iV; a)- Theories of superdense matter en- 
ter in Eq. ([T]) through the energy-momentum tensor. 



(e + P)+.g^,P, 



(2) 



which contains the equation of state (EoS), i.e. pressure as 
a function on energy density, P(e), of the stellar matter. 
The quantities and Un in Eq. ([2|) are four-vefocities, 
defined as = dx^/dr and Ui, = dx^/dT, where dr^ — 
—ds^. These velocities are the components of the macro- 
scopic velocity of the stellar matter with respect to the 
actual coordinate system that is being used. 



2.1 Non-rotating stars 
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Fig. 2. Neutron star mass versus central energy density for 
different equations of state of neutron star matter [8]. 



Non-rotating neutron stars are spherically symmetric. The 
metric of such stars thus dependes only on the radial co- 
ordinate and is given by 



ds' 



2 s> 1.2 I 2 /I J 2 
-e dt -f e dr 



'de' 



sin^e d(t)^ , (3) 



where <l> and A are the radially varying metric functions. 
From Eq. ([3]) it follows that the components of the metric 
tensor are given by 
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sin^ei , (4) 



so that the only non-vanishing Christoffel symbols are 
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rW = -r sin^e e"^"^ ,rL = - sin 9 cos ( 



Subsituting Eqs. (|4|) and (O into ^ and using 



(5) 



= 



(the semicolon denotes covariant differentiation) leads to 
the Tolman-Oppenheimer-VolkofT (TOV) equation [2j, 



dP_ 

dr 



(l-fP/e) (l + ATTr^P/m) 
1 — 1m I r 



(6) 



applicable to compact stellar configurations in general rel- 
ativistic hydrostatic equilibrium. We use units for which 
the gravitational constant and velocity of light are G = 
c = 1 so that the mass of the sun is M© = 1.47 km. The 
mass contained in a sphere of radius r is given by m = 
47r /J" r^edr. Hence the star's total mass follows as M = 
m{R), where R denotes the star's radius. The Newtonian 
limit of Eq. ^ is obtained for P/e < 1, Anr^P/m < 1, 
and 2m /r <^ 1. General relativity, thus, increases the pres- 
sure gradient inside the star which leads to smaller, more 
compact stars. The masses of neutron stars lie between 
about one and two solar masses, and their radii are around 
10 km. Thus, 2M/R - 30 - 60% for neutron stars. Solu- 
tions of the TOV equations are shown in Fig. [2l 



2.2 Rotating stars 

The structure equations of rotating neutron stars are con- 
siderably more complicated than those of non-rotating 
neutron stars [5]. These complications have their cause 
in the rotational deformation, that is, a flattening at the 
pole accompanied by a radial expansion in the equatorial 
direction, which leads to a dependence of the star's met- 
ric on the polar coordinate, 6. Secondly, rotation stabilizes 
a neutron star against gravitational collapse. A rotating 
neutron star can therefore carry more mass than a non- 
rotating star. Being more massive, however, means that 
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Fig. 3. Dragging of local inertial frames (Lense-Thirring effect) 
caused by ~ 1.4 Mq neutron stars rotating at 2 ms [9]. The 
frequency ui is defined in Eq. ((9}. 



the geometry of space-time is changed too. This renders 
the metric functions of a rotating neutron star frequency 
dependent. Finally, the general relativistic effect of the 
dragging of local inertial frames implies the occurrence of 
an additional non-diagonal term, g*'^, in the metric tensor 
(7^". This term imposes a self-consistency condition on the 
stellar structure equations, since the degree at which the 
local inertial frames are dragged along by the star is de- 
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termined by the initially unknown stellar properties like 
mass and rotational frequency. The covariant components 
of the metric tensor of a rotating compact star are thus 
given by 

gee = e^^ , g^^^e^^ , (7) 

which corresponds to a line element, ds^ = g^i^dx^dx^ , of 
the form 

ds^ = -e^^dt^ + e2^(d(/) - ujdtf + e^^^dO^ + e^^dr^ . (8) 

Here each metric function, i.e. u, ip, jj, and A, as well as the 
angular velocities of the local inertial frames, uj, depend 
on the radial coordinate r and on the polar angle 9 and, 
implicitly, on the star's angular velocity f2. Of particular 
interest is the relative angular frame dragging frequency, 
ui, defined as 

w(r, 0, i7) = i7 - w(r, 6*, J7) , (9) 

which is the angular velocity of the star, J7, relative to 
the angular velocity of a local inertial frame, w. It is this 
frequency that is of relevance when discussing the rota- 
tional flow of the fluid inside the star, since the mag- 
nitude of the centrifugal force acting on a fluid element 
is governed-in general relativity as well as in Newtonian 
gravitational theory-by the rate of rotation of the fluid 
element relative to a local inertial frame [TT] . In contrast 
to Newtonian theory, however, the inertial frames inside 
(and outside) a general relativistic fluid are not at rest 
with respect to the distant stars. Rather, the local inertial 
frames are dragged along by the rotating fluid. Depend- 
ing on the internal stellar constitution, this effect can be 
quite strong, as shown in Fig. [3] for neutron stars rotat- 
ing at 2 ms [9] . The stellar models are computed for three 
different equations of state: (1) HV, which describes neu- 
tron stars made of nucleons (N) and hyperons (H) ; Gf ° , 
which describes neutron star matter in terms of nucleons, 
hyperons, and quarks (Q); and CFL, which assumes that 
the quarks are color superconducting (CFL phase) [H]. 
For a very compact neutron star, as obtained for the CFL 



case, one sees that the local inertial frames at the star's 
center rotate at about half the star's rotational frequency, 
uj{r = 0) ~ J7/2. This value drops to about 15% for the 
local inertial frames located at the star's equator. 

Figures [4] and [5] show the influence of rotation on the 
mass-radius relationship of neutron stars. For ultrafast 
rotation at the Kepler frequency, a mass increase up to ~ 
20% is obtained, depending on the equation of state. The 
equatorial radius increases by several kilometers, while the 
polar radius get smaller by several kilometers. The ratio 
between both radii is around 2/3, except for rotation close 
to the Kepler frequency. 

2.3 Limiting rotational periods 

No simple stability criteria are known for rapidly rotat- 
ing stellar configurations in general relativity. However, 
an absolute limit on rapid rotation is set by the onset 
of mass shedding from the equator of a rotating star. The 
corresponding rotational frequency is known as the Kepler 
frequency, f^K- In classical mechanics, the expression for 
the Kepler frequency, determined by the equality between 
the centrifugal force and gravity, is readily obtained as 
J7k = \/ M/B?. In order to derive the general relativistic 
counterpart of this relation, one applies the extremal prin- 
ciple to the circular orbit of a point mass rotating at the 
star's equator. Since r = 9 = const for a point mass there, 
one has dr = d9 = 0. The line element ([S]) then reduces to 
ds"^ — {e^" — e^^ {Q — ujY) dt^ . Substituting this expression 
into J = ds, where si and S2 refer to points located at 
that particular orbit for which J becomes extremal, gives 



dt^e^^ -e'^-^ {n-ujf . (10) 

J si 

Applying the extremal condition (5J — to Eq. (fTO|l and 
noticing that V = e^^^ (Q — uj) then leads to 

i^^re^'^V^ - uj^re'^+'^V -ly^re^" = 0. (11) 

This relation constitutes a quadratic equation for the or- 
bital velocity F of a particle at the star's equator. One 
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Fig. 6. Moment of inertia of several sample neutron stars. 



thus obtains for the Kepler frequency /?k the final rela- 
tion [2] 



which is to be solved self-consistently at the equator of 
a rotating neutron star. The general relativistic Kepler 
period follows from Eq. as Pk — 27r/i?K- For typ- 
ical neutron star matter equations of state, the Kepler 
period obtained for 1.4 Mq neutron stars scatters around 
1 ms. One exception to this are strange quark matter stars. 
These are self-bound and, thus, tend to possess smaller 
radii than conventional neutron stars, which are bound by 
gravity only. Because of their smaller radii, strange stars 
can withstand mass shedding down to periods of around 
0.5 ms [I31IT4]. The CFL neutron star discussed just above 
is a quark-hybrid star and as such bound by gravity only. 
The dense quark core in the center of this neutron star 
implies a relatively large binding energy of 0.12 Mq, lead- 
ing to a rather low mass shedding period of 0.7 ms. In 
closing this section, we introduce the moment of inertia 
of a rotating neutron star described by the metric in Eq. 
([5]). For such stars the moment of inertia is given by 



27r / de 



dr 6^+'"+"+'^ 



(13) 



Figure [H] shows that the crustal fraction of the moment of 
inertia of a neutron star may be around 50% smaller if the 
star contains a very soft phase of matter like CFL quark 
matter. This may be of relevance for pulsar glitch models 
and the modeling of the post-glitch behavior of pulsars. 



3 Composition of neutron star matter 

A vast number of models for the equation of state of neu- 
tron star matter has been derived in the literature over the 
years. These models can roughly be classified as follows: 

— Thomas-Fermi based models [TSllTBl 



— Schroedinger-based models (e.g. variational approach, 
Monte Carlo techniques, hole line expansion (Brueck- 
ner theory), coupled cluster method. Green function 
method) [l71[T8l[T9l[20] 

— Relativistic field-theoretical treatments (relativistic mean 
field (RMF), Hartree-Fock (RHF), standard Brueckner- 
Hartree-Fock (RBHF), density dependent RBHF (DD- 
RBHF) [Timi[22p5P?ll^[^ 

— Nambu-Jona-Lasinio (NJL) models p7l[28ll29l[30ll3Tp2] 

— Chiral SU(3) quark mean field model [33] . 

Neutron star masses computed for some of these models 
are shown in Fig. [21 



3.1 Relativistic nuclear field-theoretical models 

Relativistic nuclear field-theoretical models [T|[^ [^T1I^^ 23 1 
[24l[25l[26] are based on Lagrangians of the form C = Cb + 
J^M + ^int + J^iept, where 



(14) 



C 



1 ^ {d.'P.d'^'P, - mjl>j) 

i—(T.5 



M 



2 ^ V2 

K—LJ^p 

A„t = VS4(^, V')V"?. - ^/L(^, ^)7mV'^'"^^ 

+i^rs{i^, i;)Ti;^s - ^h^T^A'-p')^' .(16) 

Here, Cb and Cm are the free baryonic and the free mesonic 
Lagrangians, respectively, and interactions are described 
by Cint, where pjit^ = d^Ai^^ —d^A^^^ is the field strength 
tensor of one of the vector mesons (k = u, p). In the case 
of RMF, RHF and RBHF the meson-baryon vertices Fa 
(a — (T, yo) are density-independent quantities which 
are given by expressions like Fg. = ig„ for the scalar a me- 
son, = g^jt" + {i/2){f^/2m)dx[l^n^] for uj mesons, 
etc. [5]. In the framework of the DD-RBHF scheme, the 
meson-baryon vertices Fa depend on the baryon field op- 
erators [21]. The field equations that follow from Eqs. 
([T4| - ([T6l) have the mathematical form 

{il^d^~mB)i^B{x)^ J2 M{x)rM^B{x), (17) 

M=i,K 

{d^'d^ + ml)cj{x) = E ^b{x) t ^b{x) , (18) 



plus similar equations for the other mesons [2ll24j. 



3.2 Hyperons and baryon resonances 

At the densities in the interior of neutron stars, the neu- 
tron chemical potential, /u", is likely to exceed the masses, 
modified by interactions, of S, A and possibly S hyper- 
ons [2',!34j . Hence, in addition to nucleons, neutron star 
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matter may be expected to contain significant popula- 
tions of strangeness carrying hyperons [34j . The thresh- 
olds of the lightest baryon resonances {A~ , A'^, Z\++) 
are reached for relativistic mean-field (RMF) calculations 
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Fig. 7. Composition of neutron star matter in RMF. 



at densities which correspond to unstable neutron stars. 
This is different for relativistic Brueckner-Hartree-Fock 
(RBHF) calculations where Z\'s appear rather abundantly 
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Fig. 9. Hyperon composition of a rotating neutron star in 
equatorial direction. (DD-RBHF calculation performed for 
Bonn (model 1) potential, non-rotating star mass is I.TOM©.) 
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via weak reactions like 

n ^ p + e~ + (19) 

to the lower energy state. The chemical potentials associ- 
ated with reaction (|19p in equilibrium obey the relation 

A^" = M^ + M^", (20) 

where fi"^' = since the mean free path of (anti) neutrinos 
is much smaller than the radius of neutron stars. Hence 
(anti) neutrinos do not accumulate inside neutron stars. 
This is different for hot proto-neutron stars Equa- 
tion ([20]) is a special case of the general relation 



in stable neutron stars [35], compare Figs. [7] and [51 De- 
pending on the star mass and rotational frequency, the to- 
tal hyperon population in neutron stars can be very large 
[34] . which is illustrated graphically in Figs. [9HT21 for ro- 
tating neutron stars based on equations of state computed 
in the framework of the DD-RBHF formalism. The stars 
shown in these figures have rotational frequencies rang- 
ing from zero to the mass shedding frequency, vk- Pure 
neutron matter, therefore, constitutes an excited state rel- 
ative to hyperonic matter which would quickly transform 



fi^ ^ B^li"" - q^ii'' , (21) 

which holds in any system characterized by two conserved 
charges. These are in the case of neutron star matter elec- 
tric charge, g-^, and baryon number charge, B^. Applica- 
tion of Eq. ([H]) to the A hyperon [B^ = 1, = 0), for 
instance, leads to /i"^ = /i" . Ignoring particle interactions, 
the chemical potential of a relativistic particle of type % 

is given by /i-^ = Lo{kp^) = y^m^ -I- fc|,^ , where Lu{kF^) is 

the single-particle energy of the particle and kp^ its Fermi 
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Fig. 11. Hyperon composition of a rotating neutron star 
in equatorial direction. (DD-RBHF calculation performed for 
Groningen potential, non-rotating star mass is 1.60 M©.) 
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momentum. One thus obtains 



> 



3 fm 



-1 ^ _ '^F„ 

^ n — 



37r2 



2no, (22) 



where ttia — 1116 MeV and m„ = 939 MeV was used. 
That is, if interactions among the particles are ignored, 
neutrons are replaced with yl's in neutron star matter at 
densities as low as two times the density of nuclear mat- 
ter. This result is only slightly altered by the inclusion 
of particle interactions [M]. Densities of just ~ 2no are 
easily reached in the cores of neutron stars. Neutron stars 
may thus be expected to contain considerable populations 
of yl's, Z"s and S's, as confirmed by the outcome of DD- 
RBHF calculations shown graphically above. Depending 
on the star's mass, the total hyperon population can be 
very large [34] . 
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